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SUMMARY 


This  report  covers  research  and  operation  activities  at  the 
Norwegian  Seismic  Array  (NORSAR)  for  the  period  1 July  - 
31  December  1975. 

The  different  systems  at  the  array  have  performed  satisfactorily 
in  the  period  (Chapter  ID , except  for  the  Special  Processing 
System  (SPS)  which  has  shown  a more  and  more  deteriorating 
performance,  ihis  situation  is  probably  not  going  to  improve 
with  time,  and  measures  should  therefore  be  taken  to  establish 
a backup  system  and/or  to  replace  this  problematic  unit.  Ex- 
tensive tests  have  been  run  with  a new  version  of  the  Detection 
Processor  (DP) , exchanging  messages  with  the  Communications 
and  Control  Processor  located  at  the  Seismic  Data  Analysis 
Center  (SDAC) , Alexandria,  over  the  ARPANET.  The  total  number 
of  events  reported  this  period  is  somewhat  smaller  than  previous 
periods  (Chapter  III),  which  can  be  explained  by  the  lack  of 
larger  seismic  swarms.  In  the  DP  an  automatic  'possible  explo- 
sion* and  'large  event'  warning  has  been  implemented  (Chapter 
IV)  . 


The  period  is  characterized  by  a stable  and  satisfactory 
operation  of  the  array  instrumentation  (Chapter  V) , result- 
ing in  a normal  maintenance  activity.  In  this  period  14  papers, 
1 report  and  8 new  programs  have  been  finished  (Chapter  VI) . 

Chapter  VII  summarizes  the  research  activities  in  the  period, 
and  altogether  10  different  topics  are  covered.  First  a study 
of  the  upper  mantle  beneath  Fennoscandia  and  Western  Russia 
is  presented.  Then  comes  an  investigation  of  the  3-dimensional 
structure  under  LASA  and  central  California.  The  corresponding 
investigation  for  NORSAR  was  described  in  the  report  for 
the  period  1 July  1974  - 30  June  1975.  A new  contribution 

to  the  discussion  about  the  precursors  to  P'P'  is  found  in 
Section  VII. 3,  while  Section  VII. 4 describes  an  analysis 
of  Rayleigh  waves  from  Novaya  Zemlya  explosions.  Two  studies 
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of  operational  capabilities  have  been  made,  one  for  NORSAR, 
using  three  years  of  data,  and  one  for  a set  cf  Fennoscandian 
stations,  using  six  years  of  data.  The  purpose  of  the  investi- 
gation presented  in  Section  VII. 7 is  to  estimate  the  M -nu 

s b 

relationship  for  earthquakes  and  assess  the  impact  of  network 
magnitude  bias  effects  on  the  results.  A new  earthquake- 
explosion  discriminant,  replacing  Mo  with  the  long  Deriod 
power  spectrum  estimate,  is  given  in  Section  VII. 8.  Section 
VII. 9 describes  the  last  results  on  a study  of  seismicity  in 
Fennoscandia , while  Section  VII. 10  describes  a project  for 
precise  monitoring  of  seismic  velocities. 


II. 
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OPERATION  OF  ALL  SYSTEMS 

II . 1 Detection  Processor  Operation  (DP) 

The  overall  goal  of  having  minimum  system  down  time  for  the 
Detection  Processor  System  has  in  this  period  to  a small 
extent  been  in  conflict  with  the  demands  for  testing  out 
the  new  version  of  the  system  as  the  'primary  online'  system. 
This  was  especially  felt  towards  the  end  of  the  period.  Aoart 
from  this,  hardware  problems  have  mainly  been  the  causes  of 
breaks  in  the  recording.  The  up  time  percentage  is  97.3%,  as 

compared  to  99.2%  for  the  previous  reporting  period  (July  1974 
to  June  1975) . 

Fig.  II. 1.1  and  the  accompanying  Table  II. 1.1  both  show  the 
daily  DP  down  time  In  hours  for  the  days  between  1 July  and 
31  December  1975.  The  monthly  recording  times  and  up  percentages 
are  given  in  Table  II. 1.2. 

The  most  significant  break  in  the  recording  occurred  on  3 to 
4 September,  when  a hardware  error  in  the  SPS  Multiplexer 
caused  a down  period  of  more  than  41  hours.  Also,  thermal 
problems  caused  the  SPS  to  stop  26  and  27  December  for  more 
than  16  hours,  while  a faulty  core  component:  caused  an  SPS 
stop  of  21  hours  29  and  30  December.  From  Table  II. 1.2  we 
see  that  we  in  July  had  41  stops  and  a mean  time  between 
failures  of  only  0.7  days.  The  main  cause  of  these  frequent 
stops  was  that  the  DP  did  not  receive  an  interrupt  from  the  SPS 
(Special  Processing  System)  within  the  required  time  span. 

The  115  breaks  can  be  grouped  in  the  following  cateogries: 


SPS  Problems 

• 

• 

45 

Test 

Hardware  malfunction  (air  condition, 

• 

• 

18 

CPU,  1052) 

• 

• 

11 

Power  breaks  & jumps 

• 

• 

10 

Unknown  reason 

• 

• 

7 

Disk  drive  problems 

• 

• 

6 

Maintenance  (lines  or  CE)  : 5 

Tape  drive  problems  . 4 

Software  problems  . 3 

TOD  unit  problems  . 2 

ARPANET  TIP  . 2 

Program  change  . ± 

Operator  error  . ^ 

The  'Tests'  category  consists  of  relatively  small  breaks  in 
the  recording,  when  the  new,  experimental  DP  was  taken  up  as 
the  'primary'  system.  Conclusively,  it  ir>  evident  from  this 
that  malfunctioning  or  breakdown  of  the  SPS  is  the  cause  of 
most  of  the  breaks  in  recording,  and  also  causes  the  lenoest 
breaks.  It  is  not  reasonable  to  believe  that  this  situation 
will  improve  with  time.  Measures  should  therefore  be  taken 
to  counteract  an  overall  deteriorating  performance  of  the 
system.  One  way  of  doing  this  would  be  to  create  a backup 
system  capable  of  data  recording,  to  be  taken  un  when  the 
SPS  for  some  reason  went  down. 


The  total  dovrn  time  for  this  period  was  129  hours  25  minutes. 
The  mean  time  between  failures  for  the  period  was  1.6  days, 
as  compared  with  2.8  days  for  the  last  reporting  period 
(July  1974  to  Jure  1975). 

11 *2  Event  Processor  Operation  (EP) 

Apart  from  one  hangup,  caused  by  DP  writing  garbage  into  the 
Shared  Disk  Time/Tape  file  during  a re-initialization,  the 
Event  Processor  has  performed  satisfactorily  during  this 
period.  Its  up  time  percentage  is  25.5%,  as  compared  with 
29.8%  for  the  last  reporting  period  (July  1974  to  June  1975) 
(see  Table  II. 1.2).  The  computer  time  per  processed  event  is, 
however,  almost  exactly  the  same. 
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14 

14 

UPDATE  FROM  b 1 NTFREERE 

2*2 

2 

2 

2 

IV 

UNKNOWN  REAAHN 

2*5 

4 

* 

4 

27 

UNKNUWN  RLAAIJN 

3C4 

? 0 

1 C 

20 

2 V 

CHANGE  OF  TAPE  DRIVES 

3C7 

15 

16 

16 

40 

CHANGE  OF  TAPE  DRIVES 

312 

2 

13 

2 

26 

PLOT  TASK  TROUBLE 

121 

2 

4 6 

4 

1/ 

POWER  JUMP 

.12  b 

2 4 

1 5 

24 

0 

POWER  JUMP 

12  f 

0. 

0 

0 

16 

POWER  JUMP 

127 

0 

2 5 

0 

4 3 

POWER  JUMP 

12* 

1 

10 

2. 

18 

TIP  TRHlJilLL 

117 

4 

4 0 

5 

44 

POWER  JUMP 

11/ 

5 

6 4 

6 

6 

POWER  JUMP 

11  7 

9 

29 

9 

40 

DISK  ERROR 

11* 

15 

4 6 

15 

5 7 

SO AC  TEST 

149 

20 

68 

21 

4 4 

SO  AC  TEST 

1 42 

16 

2 

16 

48 

SOAC  TEST 

142 

20 

12 

20 

21 

SUAC  TESI 

1*4 

16 

2 2 

16 

47 

SHAG  TEST 

144 

18 

4 i 

18 

45 

SOAC  TEST 

1*4 

20 

4 5 

20 

40 

SOAC  TEST 

14b 

15 

56 

16 

2 

SDAC  TFST 

14b 

17 

5 

17 

21 

SOAC  TEST 

14b 

17 

61 

17 

57 

SOAC  TEST 

14b 

IV 

63 

20 

0 

SOAC  TFST 

151 

15 

6 4 

16 

0 

SUAC  TCST 

451 

'4 

1U 

14 

6 3 

SOAC  TEST 

Vjl 

/.  .4 

15 

23 

31 

SOAC  TEST 

451 

2 1 

6 * 

24 

n 

SOAC  TESI 

452 

0 

r; 

0 

13 

SDAC  TEST 

452 

C 

2 9 

0 

4 9 

SOAC  TESI 

VI 2 

15 

65 

16 

0 

SDAC  TEST 

162 

19 

l 

19 

26 

SOAC  TEST 

460 

1C 

9 

11 

10 

SPS  THFRMAL  PROBLEMS 

400 

14 

2 4 

14 

51 

SPS  THERMAL  PROBLEMS 

4 60 

16 

44 

20 

7 \ 

SPS  THLRMAI  PR 08|  r M S 

460 

22 

l‘» 

2 4 

i ; 

SPS  THERMAL  PROBLEMS 

460 

24 

4 V 

24 

c 

SPS  THERMAL  PROBLEMS 

461 

0 

t ; 

l 

26 

SPS  THERMAL  PROBLEMS 

361 

6 

4 9 

10 

3 

SPS  THF.RMA.  PROBLEMS 

461 

to 

1 1 

12 

7 

SPS  THERMAL  PROBLEMS 

361 

16 

1 4 

24 

0 

SPS  HARDWARE  EKRHR 

464 

0 

■) 

14 

1 7 

SPS  HARDWARE  ERROR 

465 

21 

4 0 

24 

0 

NEW  YEAR  STUP 

Table  II 
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II. 3 NORSAR  Data  Processing  Center  (NDPC)  Operation 
Data_Center 

No  change  in  operational  procedures,  eqi -pment  configuration 
or  facilities  occurred  in  the  report  period. 

Maintenance  of  equipment  also  continued  as  before,  mainly 
performed  by  subcontractors  for  the  different  categories 
of  equipment.  Project  personnel  perform  maintenance  on  some 
of  the  special  equipment  not  covered  by  standard  contracts. 

Equipment  performance  in  general  was  reasonably  stable,  as 
indicated  by  a Detection  Processor  up- time  of  97.3%  for  the 
period.  Of  the  remaining  2.7%,  0.2%  is  accounted  for  by 
power  failure,  2%  by  Special  Processing  System  (SPS)  failure. 
Peripheral  equipment  have  a somewhat  higher  failure  rate, 
although  usually  with  no  serious  effect  on  essential  functions. 

The  SPS  continues  to  cause  some  concern.  It  has  had  a series 
of  - so  far  - unexplained  stops.  On  one  occasion  (Dec  30) 
a faulty  component  was  found  (stops  have  also  occurred  since)  . 

It  should  be  pointed  out  that,  because  of  the  complexity 
and  uniqueness  of  this  machine,  local  expertise  is  rather 
limited. 

The  ARPANET  Terminal  Interface  Processor  (TIP)  and  the 
360/TIP  interface  units  have  performed  satisfactorily.  Minor 
modifications  to  the  TIP  were  performed  by  Bolt  Beranek  and  Newman 
Inc.  ( BBN ) personnel  in  August  and  again  in  September. 

Communications 

The  subarray  communications  system  had  about  'normal'  outage 
compared  to  the  long-term  average.  A few  stations  (04C,  07C 
and  13C)  had  longer  outages  caused  by  broken  ground  cables. 

The  high  outage  figure  for  01C  is  caused  by  power  failure 
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(broken  ground  cable)  rather  than  communications  line  fault 
(Table  II. 3.1) . 

An  extensive  program  for  remeasuring  and  reconditioning  of 
communications  circuits  to  conform  with  CCITT  M102  specifications 
in  cooperation  with  the  Norwegian  Telegraph  Administration 
(NTA) , is  in  progress. 

The  international  circuits  (London,  SDAC)  appear  to  be  stable, 
with  only  short  and  occasional  breaks  being  observed.  No 
statistics  for  these  circuits  are  being  kept,  as  no  automatic 
registration  is  being  performed  here. 
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Table  II. 3.1 

Communications,  degraded  performance  (>20) /outages  (>200). 

Figures  in  per  cent  of  total  time.  Month*  4 or  5 week  periods  as  indicated. 


Sub- 

array 


AVG.  0.3 


AUG (4) 
>20  >200 


0.3  0.1 
0.2  0.3 
0.3  0.1 
0.4  0.2 
0.9  C.b 
0.3  0.2 
0.7  0.5 


.3  0.4 


0.4  0.6 

0.3  0.9 


0.2  0.1 


0.2  0.5 
0.2  0.5 
4.9  1.9 
0.6  0.2 
0.5  0.3 
0.2  0.2 


NOV (4)  DEC (4)  AVG.  »jYR 
>20  >200  >20  >200  >20  >200 


0.6  0.1 


2.4  0.3  0.2 

P.0  0.2  0.2 

1.2  0.1 

0.4  0.2  0.1 


0.3  0.1 

0.1  0.1 


0.2  0.3 
0.1  0.1 
2.3  14.5 
C.l  0.3 


1.3  0.2  0.3  0.3  0.6 

3.2  0.3  0.5  0.5  2.1 

4.0  0.6  0.2  7.9  1.4 

14.0  0.2  0.2  0.4  2.5 

17.7  0.1  15.4  0.9  6.9 

0.9  0.3  0.2  0.2  0.3 
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II. 4 ARPANET 

During  the  reporting  period  NORSAR  has  had  the  following 
attachments  to  the  Terminal  Interface  Processor  (TIP)  located 
in  the  NORSAR  Data  Processing  Center: 

Two  terminals, 

Tektronix  4010-1  and 
Data  Dynamics  390 
Two  Special  Host  Interfaces, 

attached  to  Standard  Host  Interfaces  0 and  3 in 
the  TIP. 


The  terminals  have  routinely  been  used  for  bulletin  and  data 
transfer,  as  well  as  for  message  exchange,  the  latter  especially 
in  connection  with  testing  and  debugging  of  the  new  seismic 
data  exchange  (see  below) . The  terminal  activity  has  been 
based  upon  our  account  at  the  OFFICE-1  host  computer,  as  well 
as  an  account  on  the  SRI-AI  host  computer.  Both  these  ARPANET 
Hosts  are  located  at  Stanford,  California. 

Throughout  the  period  extensive  tests  have  been  run  with  a new 
version  of  the  Detection  Processor  (DP)  System  exchanging 
messages  with  the  Communications  and  Control  Processor  located 
at  SDAC,  Alexandria,  over  the  ARPANET.  The  new  DP  version 
incorporates  a Network  Cont'  Program  (NCP)  as  an  extra  task, 
enabling  it  to  use  ARPANEr  as  a < ommunication  medium,  instead 
of  the  old  Trans-Atlantic  v u)  line,  which  was  based  on  com- 
munication exchange  between  two  SPS  systems.  As  one  expected, 
the  tests  helped  to  track  down  various  errors  in  the  format 
and  content  of  the  exchanged  messages.  One  major  problem 
encountered  was  the  transmission  delay  related  to  the  limited 
bandwidth  available,  especially  during  peak  hours  with  heavy 
traffic.  This  delay  caused  intolerable  buildup  of  waiting 

messages  in  both  ends  and  consequent  gaps  in  the  transmission. 
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This  problem  was,  however,  solved  by  an  IMP  software  modifica- 
tion which  increased  the  available  bandwidth  sufficiently 
for  the  time  being.  Towards  the  end  of  the  period,  the  new 
DP  system  was  run  repeatedly  as  "Primary  Online",  being  able 
to  communicate  short  period  data  requests  from  SDAC  to  the 
SPS,  as  well  as  forwarding  Event  Processor  results  to  SDAC. 

Apart  from  some  bug  preventing  the  latter  function  to  be  per- 
formed correctly,  the  exchange  works  reasonably  well  for  long 
periods  of  time.  However,  the  major  remaining  problem  is 
scarcity  of  queue  storage  in  the  DP  system,  and  how  this  is  dis- 
tributed among  the  various  queue  block  sizes  at  system  start- 
up. Presently  the  situation  is  unstable,  causing  the  system  to 
terminate  when  demands  for  queue  blocks  are  critical  but 
cannot  be  fulfilled.  The  situation  might  be  improved  by  try- 
ing to  make  more  queue  storage  available,  as  well  as  preventing 
the  system  from  terminating  when  the  consequences  of  this 
does  not  influence  the  overall  performance  seriously  (e.g., 
some  non-critical  message  not  printed) . 
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III.  ARRAY  PERFORMANCE 

Some  basic  statistics  on  the  EP  operation  are  given  in 
Table  III.l,  which  shows  the  analyst  decisions  for  all  the 
DP  detections  processed  by  EP  during  the  reporting  period. 

The  percentages  are  quite  close  to  those  from  the  previous 
reporting  periods,  although  the  total  number  is  somewhat 
smaller,  which  can  be  explained  by  the  lack  of  larger  seis- 
micity swarms.  This  is  seen  from  Fig.  III.l,  where  the  number 
of  events  on  a daily  basis  is  shown.  On  a monthly  basis, 
the  results  are  given  in  Table  III. 2.  * 


With  respect  to  detection  and  location  capabilities,  we 
refer  to  Section  VII. 5. 
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Table  III.l 

Analyst  decisions  for  detections  processed  by  EP  during  the  time 
period  July-Decerbtfr  1975 


Analyst 

Classification 

Number  of 
Processings 

Percentage 

Accepted  as  events 

3214 

49.9 

Rejected  as  being 

~ Poor  SNR  or  noise 

1258 

19.6 

- Local  events 

836 

13.0 

- Double  processings 

695 

10.8 

- Communication  errors 

433 

6.7 

Sum  Processed 

~ I 

6436 

100.00 

Table  III. 2 


Number  of 

teleseismic  and  core  phase  events 
period  July-Deceraber  1975 

reported  during  the  time 

Month 

Teleseismic 

Core 

Sum 

Jul  75 

450 

221 

671 

Aug 

369 

158 

527 

Sep 

312 

83 

395 

Oct 

346 

100 

446 

Nov 

380 

84 

464 

Dec 

314 

74 

388 

190  200  210  220  230  240  250  260  270  260  290  300  310  320  330  340  350  360 


DRY  OF  YEAR 


Fig.  III.l 

4 


Number  of  events  as  a function  of  day  of  year 
July  - December  1975. 
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IV  IMPROVEMENTS  AND  MODIFICATIONS 

IV . 1 Detection  Processor 

Froiu  the  end  of  October  the  Detection  Processor  will  check 
the  reduced  Selected  Surveillance  detections  to  see  if  a 
detection  is  from  one  of  a set  of  beam  directions  (2,  55, 

117  or  126)  all  related  to  possible  explosion  test  sites. 

If  a detection  is  on  one  of  these  beams,  and  has  a signal- 
to-noise  ratio  above  a certain  threshold,  a 'Possible 
Explosion'  warning  message  will  be  printed  on  the  1052 
Console  Typewriter  as  well  as  on  the  1403  Printer,  to  alert 
the  operator.  Also,  if  the  MSTA-value  of  a Selected  Surveillance 
detection  exceeds  a certain  threshold,  a 'Large  Event'  warning 
message  will  be  printed  on  the  1052  Typewriter  and  the  1403 
Printer.  Detailed  instructions  have  been  worked  out  and  imple- 
mented for  the  operators  when  these  messages  occur.  The  basic 
procedure  is  that  the  Strip  Chart  Recorder  is  started  with  on- 
line long  period  data  as  soon  as  the  message  occurs,  the  event 
is  immediately  processed  on  Fast  EP,  and  a seismologist  is 
called  upon  in  accordance  with  a certain  system  of  priorities. 


IV. 2 Event  Processor 

No  modifications  were  performed  in  the  EP  system  during 
this  period. 


% > 
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IV . 3 Array  Instrumentation 

The  status  of  a number  of 

Depression  of  noise 
in  SLEM  discrete 
inputs* 

Too  low  surge  rating 
of  BE  protection  cards* 

False  triggering  of 
CTV  water  monitor 

Trends  towards  nega- 
tive DC  offset  in  the 
SP/LTA 


improvement  projects  is  as  follows: 

The  material  needed  for  modification 
has  been  procured  and  a few  modules 
are  ready  for  field  installation. 

Thirty-one  cards  are  modified  and 
ready  for  field  installation. 

The  "turn-on  point"  has  been  ad- 
justed on  all  subarrays. 

Modification  (refer  Larsen  and 
Nilsen,  1974)  is  accomplished  on 
the  last  five  subarrays  left  over 
from  last  reporting  period. 

Pettersen,  1975. 


* Refer  Larsen,  Falch  and 

At  14C  the  unit  installed  on  channel  05  29  April  1975  dividing 
the  net  frequency  18  times  to  2.78  Hz  with  signal  amplitude 
of  9.2  volts  peak-to-peak,  is  still  used  for  investigation 
purposes . 


As  of  19  September  1975  14C03  and  the  14C  LP  channels  were 
attenuated  -30  dB  in  order  to  receive  unsaturated  seismic 
signals  from  events  up  to  Mfa  7.0  at  a 30  degree  distance. 
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V MAINTENANCE  ACTIVITY 

This  section  includes  a review  of  the  maintenance  accomplished 
at  the  subarrays  by  the  field  technicians  as  a result  of  the 
remote  array  monitoring  and  visual  inspections.  There  are  no 
changes  in  the  monitoring  schedule  in  the  period. 

Maintenance  Visits 

Fig.  V . 1 shows  the  number  of  visits  to  the  different  subarrays 
in  the  period.  Excluding  visits  caused  by  troubles  in  the 
communication  system  and  measurements  to  verify  CCITT  M102 
line  specifications,  the  subarrays  have  on  the  average  been 
visited  4.6  times. 


No.  of 
Visits 


01B  03B  05B  07B  02C  04C  06c  OGc  10C  12C  14C 


Fig.  V.  1 Number  of  maintenance  visits  to  the  NORSAR  subarrays 
1 July  1975  - 31  December  1975. 


r » mui 


! 


■ 
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Preventive  Maintenance  Projects 

The  preventive  maintenance  work  in  the  array  is  descr 
in  Table  V.l. 


Table  V.l 


Unit 

Action 

No.  of 
ACC  . 

Channels 

Remain. 

Com 

SP 

Seis. 

Adjustment  of 
damping 

2 

12C04 , 0 

LTA 

Adjustment  of  SP  DC 
offset  to  pos.bias 

49 

Modification  of  SP 
DC  offset  adjust- 
ment range  (R12)* 

30 

04B, 01C 
13C , 14 

RA-5 

Replacement  of 
RA-5 ' s with  new 
batteries  to  pre- 
vent power  decay 

11 

12C01-0 

14C 

CTV 

Cleaning  of  CTV 

15 

(sub- 

arrays) 

01A, 01B 
05B,  01( 
04C,  06( 
08C , 09( 
11C,  13C 

* Refer  Larsen  and  Nilsen,  1974. 

Unit 


Characteristic 


Se.is-  Damping 

mometer  Nat.  Freq. 

Sensitivity 

RCD 

Seis-  Gain 

mometer  Distortion 
Ampli-  Cal. amp. inop. 

fier  Balance 

(RA-5)  Taper  pin  block 

LTA  Ch.  gain 

Filter  discr. 

DCO 

CMR 

K1  relay  fault 

BE  card 

WHV 

Prot. 

Card 

SLEM 
BB  gen. 

RSA/ADC 

EP'J 
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Malfunction  of  Rectifiers , Power  Loss,  Cable  Breakages 

Malfunction  of  the  rectifiers  and  power  loss  requiring  action 
of  the  field  technicians  or  local  power  company  are  reported 
in  Table  V.3. 


Table  V.3 

Faults  disclosed  in  subarray  rectifiers  and  power  loss 


Sub- 

array 

Fault 

— 
Period  of 
inoperation 

Comments 

01B 

Rectifier  trafo  M2  burned 

— | 

14-15  July 

02B 

Main  AC  power  break 

7-9  Dec 

02C 

Rectifier  continuous  in 
high  charge 

$ 

Timer  defective 

03C 

Rectifier  failed  to  go 
into  high  charge 

Timer  faulty 

Sixteen  cable  breakages  have  been  repaired  in  the  period  re- 
quiring 53  days'  work  of  the  field  technicians.  The  local 
power  company  completed  the  repair  of  the  power  cable  to  01C 
and  1 KV  installations  (previously  10  KV)  29  August. 

Conclusion 

The  period  is  characterized  by  a stable  and  satisfactory  opera- 
tion of  the  array  instrumentation,  resulting  in  a normal 
maintenance  activity.  The  number  of  replaced  LP  RCDs  was 
eleven  during  the  first  three  months  of  the  period,  in  the 
last  three  months  none  have  been  replaced,  thus  indicating  a 
more  stable  trend. 


Modification  of  SP  DC  offset  adjustment  range  is  now  completed. 
The  array  average  DC  offset  of  the  SP  channels  was  plus  one 
millivolt  at  the  end  of  the  period. 

Also  the  preventive  maintenance  program  for  the  WHV's  and 
RA-5's  is  completed.  This  program  is  planned  reinitiated  in 
summer  1977. 
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A number  of  BE  lightning  protection  cards  have  been  modified 
for  improved  surge  rating  and  are  ready  for  field  installation. 
This  program  is  planned  completed  on  most  of  the  subarrays  before 
the  next  lightning  season  in  summer  1976. 
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ABBREVIATIONS 


AC 

- 

Alternate  Current 

ADC 

- 

Analog-to-Digital  Converter 

BB 

- 

Broad  Band 

BE  card 

- 

Lightning  Protection  Card 

CMR 

- 

Common  Mode  Rejection 

CTV 

- 

Central  Terminal  Vault 

DC 

- 

Direct  Current 

DC0 

- 

DC  of  Etft 

EPU 

- 

External  Power  Unit 

LP 

- 

Long  Period 

LTA 

- 

Line  Terminating  Amplifier 

RA-5 

- 

SP  Seismometer  Amplifier 

RCD 

- 

Remote  Centering  Device 

RSA 

- 

Range  Switching  Amplifier 

SLEM 

- 

Seismic  Short  and  Long  Period  Electronics  Module 

SP 

- 

Short  Period 

WHV 

- 

Nell  Head  Vault 
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VI . 2 Program  Documentation 

The  following  programs  have  been  developed  at  NORSAR  during 
the  reporting  period: 

A new  improved  magnetic  tape  read  routine  that  will  allow 
the  programmer  to  build  decisions  into  his  program  about 
how  to  deal  with  data  checks,  etc.  (N/PD-82,  READTA) 

An  improved  version  of  t ie  HYPERMAP  map  program,  allowing 
more  detailed  mapping  in  11  different  projections,  for 
all  or  parts  of  the  world  (N/PD-83,  SUPERMAP) 

A program  that  calculates  long  periodic  beam  power  as 
functior  of  azimuth,  slowness  and  time  (N/PD-84,  LPBEAM) 

- A program  that  calculates  a long  periodic  array  beam, 

as  well  as  uncorrected  and  response  corrected  power  spectra 
for  individual  channels  and  the  beam  (N/PD-85,  LPBMSP) 

A program  performing  a frequency-time  analysis  of  a 
given  time  series,  by  means  of  multiple  narrow  band 
filtering  in  t.ie  frequency  domain  (N/PD-86,  FTAN) 

A program  for  plotting  individual  short  period  traces 
from  the  HR-tape,  with  trace  spacing  proportional  to 
the  projection  of  the  individual  instrument  locations 
onto  the  arrival  azimuth  (N/PD-87,  SPSECT) 

A general  purpose  subroutine  for  plotting  contours  of 
a two  dimensional  scalar  field,  represented  as  a two- 
dimensional  grid  of  N by  M observations  (N/PD-88,  CPLOT) 

A routine  giving  back  the  number  of  characters  plotted 
when  using  the  NUMBER  plotting  subroutine  (N/PD-89,  NUMCAR) . 

In  addition,  various  other  programs  have  been  either  developed 
at  NORSAR,  or  received  from  other  sources  and  modified  for  use 
here. 
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VII  SUMMARY  OF  SPECIAL  AND  TECHNICAL  REPORTS/ 

PAPERS  PREPARED 

VII. 1 P-wavc  Velocities  in  the  Upper  Mantle  beneath  Fennoscandia 
and  Western  Russia 

A detailed  and  extensive  record  section  which  has  been  con- 
structed from  recordings  at  NORSAR  of  presumed  explosions 
in  continental  Russia  exhibits  two  distinct  (T,A)  triplications. 
A particularly  noteworthy  feature  of  the  (T,  A)  curves  defined 
by  the  record  section  is  the  pronounced  extension  of  the 
first-arrival  branch  for  A<(21°  as  a secondary  arrival  to  a 
distance  A~33°.  The  reliable  identification  of  these  upper 
mantle  travel-time  branches  has  been  possible  because  of  the 
extremely  dense,  areal  sampling  of  the  NORSAR  configuration. 

The  details  of  the  P-wave  multibranching  offer  powerful  con- 
straints on  the  nature  of  the  upper  mantle  transition  zone 
beneath  the  Baltic  Shield  and  Russian  platform.  Notwithstanding 
the  questions  of  uniqueness  and  baseline  uncertainty,  the 
following  features  appear  to  be  demanded  by  the  observational 
data: 

1.  A relatively  small  positive  velocity  gradient  below  the 
Moho  to  about  420  km  depth 

2.  A relatively  rapid  increase  in  velocity  by  about  7%  near 
420  km  depth 

3.  A relatively  rapid  increase  in  velocity  by  about  4% 
near  690  km  depth 

4.  A relatively  large  positive  gradient  between  the 
'discontinuities'  near  420  and  690  km  depth. 

A simple  representative  velocity  model  KCA  has  been  derived 
which  accounts  adequately  for  all  the  principal  features 
of  the  observational  data.  In  the  interests  of  simplicity, 
the  model  KCA  includes  first  order  discontinuities  in  prefer- 
ence to  the  depth-distributed  discontinuities  of  comparable 
effect.  The  KCA  model  and  the  associated  reduced  (T, A)  curves 
are  shown  in  Fig.  VII. 1.1.  It  is  apparent  in  this  figure 
that  the  differences  in  the  model  and  inferred  observational 


A (km) 


Fig.  VII. 1.1  p velocity  model  KCA,  and  corresponding  reduced  travel 
time-distance  curves  (full  lines).  Differences  in  the 
observed  and  model  reduced  times  are  shown  by  dashed 
lines.  The  spacing  of  the  +'s  on  the  model  (T,A)  curves 
is  proportional  to  d2T/dA2  and  is  indicative  of  relative 
amplitudes  (considering  geometrical  spreading  only) . 
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(T, A)  curves  are  minor.  In  addition,  the  KCA  model  is  broadly 
consistent  with  observed  relative  amplitudes  of  the  various 
branches,  although  in  this  respect  the  uncertainties  are,  of 
course,  relatively  large  because  of  possible  departures  from 
perfect  elasticity  and  other  effects. 

The  discussion  by  Simpson  et  al  (1974)  on  the  implication s 
of  certain  velocity  models  on  the  composition  and  state  of 
the  upper  mantle  Is  fully  applicable  with  regard  to  the  KCA 
model . 

D.W.  King 

G.  Calganile  (Bari,  Italy) 
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Determination  of  the  3-dimensional  Seismic  Structure  of 
the  T.ithosphere  under  Montana  LAS  A and  under  the  USGS 


Central  California  Seismic  Array 

The  novel  3-dimensional  earth-modelling  technique  presented 
by  Aki,  Christof fersson  and  Husebye  (1976)  has  been  used  for 
detailed  investigations  of  the  seismic  structures  beneath 
the  LASA  and  Central  California  arrays.  The  characteristic 
feature  of  the  above  technique  is  flexibility  in  modelling 
the  lithosphere.  The  starting  point  ir.  the  layered  medium  of 
classic  seismology  but  each  layer  is  divided  into  many 
blocks  and  assigned  a parameter  to  each  block  which  describes 
the  velocity  perturbation  from  the  average  for  the  layer. 

The  data  used  are  teleseismic  travel  time  residuals  observed 
at  the  LASA  array  in  Montana  and  USGS  Central  California 
array.  By  isolating  various  sources  of  errors  and  biases,  we 
arrive  at  a system  of  equations  to  determine  the  model 
parameters,  i.e.,  the  velocity  perturbations  within  the 
individual  blocks.  The  final  solutions  were  obtained  by  the 
use  of  the  generalized  inverse  and  stochastic  methods  in- 
cluding analysis  of  resolution  and  errors  of  the  estimated 
parameters. 


The  lithosphere  beneath  LASA 

The  seismic  velocity  structure  beneath  LASA  has  intrigued 
seismologists  ever  since  this  large  aperture  array  was 
established  in  Montana  in  1965.  A commonly  adapted  model 
here  is  one  where  material  properties  within  layers  are 
homogeneous  but  in  which  the  crust  thickens  sharply  from 
approx.  50  km  to  approx.  60  km  towards  northwest  near  the 
center  of  the  array.  Our  results,  (Aki  et  al,  1976)  a 3- 
dimensional  seismic  image  of  the  lithosphere,  are  displayed 
in  Figs.  VII. 2.1  and  VII. 2. 2 in  the  form  of  two  cross- 

sections  for  both  the  stochastic  inverse  and  generalized 
inverse  solutions  respectively. 
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The  most  conspicuous  feature  of  the  seismic  structure  under 
LASA  is  a low  velocity  anomaly  in  the  central  and  NE  part 
of  the  array  siting  area  with  a N60°E  trend  and  persisting 
from  the  upper  crust  to  depths  greater  than  100  km.  We  are 
tempted  to  associate  this  low  velocity  zone  with  wrench 
faulting  and  shearing  during  the  Nevadan  and  Laramide  oro- 
genies. The  sharp  change  in  crustal  thickness  proposed  in 
earlier  studies  is  likely  to  be  significantly  overestimated 
because  the  upper  crust  and  the  lower  lithosphere  appear  to 
share  the  same  anomaly  pattern  as  the  proposed  Moho  topography. 

B.  The  lithosphere  beneath  Central  California  - San  Andreas  fault 

We  have  also  studied  the  seismic  structure  of  the  lithosphere 
beneath  central  California  using  data  from  the  26  U.S. 
Geological  Survey  stations  which  are  shown  in  Fig.  VII. 2. 3 
together  with  the  major  tectonic  features  of  the  area 
(Husebye  et  al,  1976). 

The  San  Andreas  fault  system  is  a part  of  the  present  boundary 
between  the  North  American  and  Pacific  plates  and  expresses 
their  relative  motions  as  a transform  fault.  Between  29  and 
21  m.y.  the  Farallon  plate  broke  up  between  the  Mendocino 
and  Murray  fracture  zones  and  the  small  pieces  under- thrust 
the  west  coast  of  North  America  (Fig.  VII. 2. 4).  The  Pacific 
Ocean  floor,  which  came  in  contact  with  the  North  American 
plate  10-20  m.y.  ago,  was  very  young  so  that  the  lithosphere 
was  weak  and  susceptible  to  deformation.  As  the  eastern 
margin  of  the  Pacific  lithosphere  cooled,  thickened  and  grew 
stronger,  it  probably  became  strong  enough  to  bite  off  a 
piece  of  California  from  North  America. 

With  this  background  information,  let  us  look  at  our  results 
for  discussion.  We  start  with  Layer  1 or  the  crust  with  an 
assumed  thickness  of  25  km.  The  trends  in  the  velocity 
anomalies  in  Layer  1 (Fig.  VII. 2. 5)  run  nearly  parallel  to 
the  direction  of  th«;»  major  fault  lines  in  the  area;  high 
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velocity  to  the  west  of  the  fault  zone  and  low  velocity  to 
the  east.  The  maximum  velocity  contrast  between  the  two 
sides  for  both  solutions  is  around  10%,  which  is  generally 
consistent  with  seismic  profiling  results  on  velocity  dif- 
ferences between  the  granitic  rocks  of  the  Gabilan  Range  and 
the  Franciscan  rocks  of  the  Diablo  Range.  We  also  notice 
that  the  velocity  variation  across  the  fault  becomes  gradual 
where  the  San  Andreas  fault  system  branches  and  sharpens 
where  the  faults  merge  to  a narrow  zone  near  Hollister. 

The  slowness  anomalies  in  Layer  2 (25-50  km)  directly  beneath 
Iloho  are  shown  in  Fig.  VII. 2. 6 for  the  generalized  inverse 
solution  which  is  different  from  the  suchastic  inverse  solu- 
tion. Because  of  the  difference  in  the  two  solutions,  we 
shall  refrain  from  any  serious  attempt  to  interpret  Layer  2 
anomalies,  but  if  the  general  trend  in  the  generalized 
inverse  solution  is  correct,  we  find  high  velocity  on  the 
east  and  low  velocity  on  the  west  side  of  the  fault,  a 
reversed  picture  as  compared  to  Layer  1 . 

The  slowness  anomalies  for  Layer  3 (50-75  km)  are  given  in 
Fig.  VII. 2. 7 showing  a definite  trend  parallel  to  the  fault 
zone.  Just  like  Layer  1,  Layer  3 has  high  velocity  to  the 
west  of  the  fault  zone,  and  low  velocity  to  the  east.  The 
remarkable  correlation  of  velocity  anomalies  in  Layers  1 
through  3 with  the  San  Andreas  suggests  that  the  fault  zone 
penetrates  to  a depth  of  at  least  75  km. 

The  slowness  anomaly  patterns  for  Layer  4 (75-100  km)  and 
Layer  5 (100-125  km)  no  longer  correlate  clearly  with  the 
fault  zone,  as  shown  in  Figs.  VII. 2. 8 and  VII. 2. 9.  The 
similarity  of  their  patterns  to  Layer  3 is  only  retained  in 
the  northern  part  of  the  area*  where  the  velocity  is  still 
high  to  the  west  and  ?ow  to  the  east.  This  suggests  that 
these  layers  may  not  belong  to  the  lithosphere,  which  is 


consistent  with  the  average  thickness  of  the  lithosphere  of 
around  80  km  in  the  western  United  States  as  obtained  b\  Lee 
and  Solomon  (1975)  from  surface  wave  attenuation  data. 

One  of  the  most  interesting  aspects  of  the  3-dimensional 
velocity  anomalies  is  the  possibility  that  the  shape  of  the 
anomaly  may  indicate  the  past  or  present  mode  of  stress  in 
the  lithosphere.  We  see  an  indication  of  a sheared  structure 
in  the  velocity  anomaly  pattern  along  the  NW-SE  profile  in 
the  Pacific  plate  west  of  the  San  Andreas  fault.  It  is 
expressed  in  a north-westward  migration  of  a high  velocity 
anomaly  along  the  Pacific  coast  with  increasing  dep«_h.  The 
pattern  is  clearer  in  the  stochastic  inverse  than  in  the 
generalized  inverse  solution.  The  direction  of  shear  implied 
by  this  pattern  may  support  the  hypothesis  that  the  convec- 
tion current  in  the  asther.osphere  is  driving  the  plate 
motion.  We  cannot,  however,  forward  a definitive  answer  to 
this  problem  because  the  profile  runs  through  the  peripheral 
region  where  the  resolution  is  poor  in  the  stochastic  in- 
verse solution  and  the  random  error  effect  is  strong  in  the 
generalized  inverse  solution.  Additional  data  currently 
being  collected  will  hopefully  give  more  definitive  results. 

Synopsi  prepared  by  E.S.  Husebye. 

K.  Aki  (M.I.T.,  Cambridge,  Mass.) 

A.  Christoff ersson  (Uppsala) 

E.S.  Husebye 
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CENTRAL  CALIFORNIA  NETWORK 


Fig.  VII. 2. 3 U.S.  Geological  Survey  Central  California  Seismograph 
stations  used  in  the  present  analysis. 


SAN  ANDREAS 


Fig.  VII. 2.4 


A schematic  cross-section  of  hypothetical  plate  boundari 
reproduced  from  Solomon  and  Butler  (1974) . 
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Fig.  VII. 2. 5 


Generalized  inverse  solution  for  Layer  1.  The  numbers 
show  the  fractional  slowness  perturbation  in  per  cent  of 
the  average  layer  slowness.  The  heavy  shaded  areas 
correspond  to  the  magnitude  of  solution  greater  than 
twice  the  standard  error.  The  letters  L and  H refer  to 
lev  and  high  velocity  anomalies  respectively.  Blocks 
not  properly  resolved  are  encircled.  The  network  center 
is  marked  by  a larger  solid  circle. 
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VII. 3 Origins  of  Precursors  to  P'P' 

Observed  precursors  to  the  seismic  phase  P'P*  consist  of 
trains  of  waves  arriving  up  to  about  150  seconds  before  the 
main  P'P'  phases.  Until  recently  these  precuraors  have 
usually  been  interpreted  as  'P'dP'  waves'  resulting  from 
reflections  of  PKP  waves  from  the  undersides  oi  postulated 
horizons  at  various  depths  d below  the  Earth's  surface  in 
regions  where  the  associated  P'P'  reflections  occur.  The 
earliest  such  P'dP'  arrivals  consistently  observed  correspond 
to  a reflecting  horizon  at  a depth  of  about  650  km  and  so 
these  precursors  are  designated  P'650P'.  In  this  particular 
case  the  evidence  supporting  the  P'dP'  interpretation,  in- 
cluding evidence  obtained  in  the  present  study,  is  extremely 
strong.  For  many  of  the  precursor  arrivals  following  P'650P', 
however,  the  P'dP'  interpretation  is  questionable. 

Seismic  scattering  in  the  lowest  200  km  of  the  mantle  (region 
D")  has  previously  been  shown  to  account  for  observed  pre- 
cursors to  PKIKP  (Haddon  and  Cleary,  1974,*  King  et  al,  1973,* 
Husebye  et  al,  1976)  while  scattering  in  the  crust  and  upper 
mantle  has  been  shown  to  account  for  observed  precursors  to 
PP  (King  et  al,  1975)  . Since  P'P'  waves  ar**  imply  PKP  waves 
which  have  undergone  reflection  at  the  ear  .r  s free  surface, 
scattering  effects  similar  to  those  associated  with  PKP  and 
PP  are  also  to  be  expected  for  P'P'. 

The  consequences  of  postulating  irregularities  in  layers  at 
the  top  and  bottom  of:  the  mantle  are  examined  theoretically 
in  respect  of  scattering  associated  with  the  usual  P'P' 
phases.  The  theoretical  results  obtained  are  then  compared 
with  detailed  observational  evidence  on  P'P'  precursor  wave- 
trains  starting  up  to  50  seconds  prior  to  the  main  P'P' 
phases  recorded  at  the  NORSAR  array. 
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The  analysis  procedure  adopted  for  investigating  P'P'  pre- 
cursor wavetrains,  called  BEAMAN  analysis,  has  been  described 
in  an  earlier  report  (see  also  King  et  al,  1975).  Some 
typical  examples  of  the  power  diagrams  obtained  by  the  BEAMAN 
procedure  are  shown  in  Fig.  VII. 3.1.  The  numbers  1,  2,  3,... 
on  the  contours  in  each  figure  specify  the  power  in  dB  below 
the  reference  level  for  the  particular  time  interval  to  which 
the  diagram  belongs.  Note  the  good  agreement  between  the 

locations  of  the  observed  maximum  power  peaks  and  the  theoretical 
results. 

The  hypothesis  of  seismic  scattering  by  small-scale  random 
inhomogeneities  in  the  uppermost  few  hundred  km  of  the  mantle 
(in  and  near  the  crust)  and  in  the  lowermost  few  hundred 
km  of  the  mantle  (region  D")  is  shown  to  be  consistent  with 
detailed  observational  data  on  precursors  to  P'P'  recorded 
at  the  NORSAR  seismic  array  for  arrival  times  0 to  50  seconds 
before  P'P' (BC) . The  data  are  generally  inconsistent  with 
theoretical  results  for  the  P'dP'  interpretation  in  this 
time  interval.  The  NORSAR  data  appear  to  be  quite  character- 
istic of  similar  but  less  detailed  data  on  precursors  to  P'P' 
obtained  by  other  authors,  previously  interpreted  as  indicating 
the  presence  of  a number  of  sharp  reflecting  discontinuities 
in  the  uppermost  200  km  of  the  mantle  below  the  crust.  The 
scattering  interpretation  removes  any  need  for  postulating 
such  discontinuities  in  this  region. 

R.A.W.  Haddon 

E.S.  Husebye 

D.W.  King 
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VII. 4 An  Analysis  of  Rayleigh  Waves  from  Novaya-Zemlya  Explosions 

The  Russian  test  sites  on  Novaya-Zemlya  are  located  in  the 
middle  of  the  island  (azimuth  * 35°  at  NORSAR)  and  on  the 
southern  end  of  it  (azimuth  * 42°) . NORSAR  has  recordings 
with  strong  Rayleigh  waves  of  6 events  from  the  northern 
test  site  and  3 from  the  southern.  For  all  these  events 
a beam  made  of  the  22  long  period  vertical  instrument 
recordings  has  been  analyzed  for  dispersion  characteristics. 

The  analysis  is  done  by  means  of  multiple  narrow  band  filtering 
in  the  frequency  domain.  The  main  output  of  the  program  is 
a map  of  power  distribution  in  the  frequency  time  domain. 

Fig.  VII. 4.1  shows  such  a map  for  an  event  in  the  northern  test 
site.  The  thick  line  on  the  figure  connects  points  oi  max in  urn 
energy.  The  most  notable  effect  on  the  figure  is  the  low  group 
velocity  observed,  say  around  2.5  km/ser  for  waves  with 
period  less  than  20  seconds.  On  Fig.  VI I. 4. 2 is  then  shown  the 
same  type  of  figure  for  an  event  from  the  southern  test  site, 
and  a much  more  complicated  pattern  arises.  First  it  is 
noted  that  for  waves  with  period  less  than  20  seconds  the 
main  (first)  arrival  now  has  a group  velocity  ol:  approximately 
2.9  km/sec.  Secondly  it  is  noted  that  two  late"  rivals 
may  be  observed  with  group  velocit-  2.5  km/sec  i 2.2  km/sec 
respectively.  Checking  the  azimui.  ;or  the  different  arrivals 
(this  has  been  done  both  with  an  ordinary  time  domain  beam 
power  program  and  with  a high  resolution  program) , it  is 
found  that  the  event  on  Fig.  VII. 4.1  has  an  azimuth  close 
to  the  correct  one.  For  the  data  on  Fig.  VII. 4. 2,  the  first  and 
third  arrivals  have  an  azimuth  close  to  the  correct  one  for 
this  event,  while  the  second  arrival  has  an  azimuth  which  is 
roughly  the  same  as  that  for  the  event  on  Fig.  VII. 4.1.  These 
preliminary  studies  indicate  that  the  extremely  low  group 
velocities  observed  for  events  from  the  northern  test  site  are 
caused  by  thick  sediment  layers  in  the  Barents  Sea  on  this 
path,  while  the  southern  path  has  somewhat  thinner  sediment 
layers.  The  multipathing  (arrival  2 and  3 on  Fig.  VII. 4. 2)  is 
believed  to  be  caused  by  discontinuities  within  the  Barents 
Sea,  the  details  of  these  will  be  studied  further. 

A.L.  Levshin  (Moscow) 

K.A.  Berteussen 


VII. 5 NORSAR  Operational  Capabilities 

Using  the  same  regionalization  and  estimation  procedures  as 
Bungum  and  Husebye  (1974)  n evaluation  of  the  NORSAR  Event 
Processor  performance  fo.  : period  April  1972  until  March 
1975  has  been  undertaken. 

Detectability  Thresholds 

Table  VII. 5.1  shows  the  estimated  50  and  90  per  cent  cumulative 
detection  thresholds  for  the  period  April  1972  - March  1973, 
April  1973  - March  1974,  April  1974  - March  1975  and  April 
1973  - March  1975.  There  is  a small  tendency  for  higher 
thresholds  for  the  last  two  years,  although  for  regions  6, 

7 and  8 the  thresholds  for  the  period  1974-75  are  lower  than 
those  for  1972-73.  The  main  impression  is  that  there  is  a 
remarkable  stability  in  the  performance  over  the  three  years. 
Using  only  one  year  intervals,  there  are  some  regions  where 
there  are  too  few  events,  or  the  distribution  is  such  that 
the  method  applied  is  not  usable.  With  three  years  of  data 
a relatively  reliable  estimate  of  the  detection  threshold 
for  all  the  regions  is,  however,  achieved. 

Location  Estimation 

The  difference  in  epicenter  solutions  has  been  calculated 
for  the  events  reported  jointly  by  USGS  and  NORSAR  in  the 
period  January  1973  until  March  1975.  Fig.  VII. 5.1  presents  the 
results  for  region  14  (distance  limits  30-90  decrees  from 
NORSAR)  in  increments  of  50  km  location  difference.  Fig.  VII. 5. 2 
shows  the  same  data  for  the  interval  February  to  November 
1972.  It  is  seen  that  the  distribution  in  Fig.vil  .5.1  has  a 
much  shorter  tail,  indicating  that  most  of  the  large  location 
errors  have  been  removed.  Since  the  data  is  so  skew  the 
50  per  cent  (median)  and  the  90  per  cent  level  of  location 

differences  is  used  as  characterizing  parameters.  In  Table 
VII. 5. 2 the  values  observed  for  the  period  February  to  November 
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1972  (by  Bungum  and  Husebye,  1974)  is  listed  together  with 
the  values  observed  for  the  period  January  1973  until  March 
1975.  For  most  of  the  regions  there  is  a considerable 
reduction  in  location  differences  especially  on  the  90  per 
cent  level.  This  is  partly  due  to  the  implementation  of  new 
region  corrections  at  NORSAR  in  November  1972,  and  partly 
due  to  the  more  experienced  analysts. 

A more  thorough  discussion  of  both  detection  and  location 
performance  will  be  published  in  a forthcoming  report. 

K.A.  Berteussen 
H.  Bungum 
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VI1-6  Event  Detectability  of  Seismograph  Stations  in  Fennoscandia 

Based  on  six  years  of  ISC  data  (1964—69)  the  operational  event 
d*  4 ection  capabilities  have  been  analyzed  for  those  stations 
in  Fennoscandia  that  reported  regularly  to  the  ISC  during 
this  time  period  (Pirhonen  et  al,  1976).  It  has  been  demon- 
strated  that  the  incremental  P-wave  detection  probability  of 
a station  as  a function  of  ISC  m^  can  be  very  accurately 
modelled  by  a cumulative  Gaussian  distribution  function. 

This  has  been  found  to  hold  true  when  analyzing  large  as 
well  as  small  regions.  By  maximum  likelihood  estimation 
based  on  the  Gaussian  model  (Ringdal,  1975)  we  have  established 
that  operational  detection  thresholds,  both  at  the  50  and  90 
per  cent  levels,  exhibit  large  variations  between  individual 
Nordic  stations.  For  the  entire  teleseismic  region,  the  50 
per  cent  thresholds  range  from  mb=4.5  to  5.5,#  the  best 
-stations  being  located  in  Finland.  This  is  consistent  with 
the  expected  increase  in  seismic  noise  level  with  decreasing 
distance  from  the  North  Atlantic  Ocean.  Fig.  VII. 6.1  shows  the 
50  per  cent  thresholds  of  all  stations  in  the  teleseismic 
range,  while  Fig.  VII. 6. 2 shows  the  estimated  detection  curves. 

Studies  of  the  noise  power  spectra  of  four  Fennoscandian 
stations  over  a one-year  period  show  as  a common  feature  a 
very  sharp  fall-off  with  increasing  frequency  in  the  band  of 
interest  for  P-wave  detection.  From  1 to  2 Hz  the  decrease 
in  average  noise  power  is  15-20  dB,  while  the  standard 
deviation  at  each  frequency  in  this  band  is  much  lower, 
typically  4-6  dB. 

An  attempt  to  predict  relative  detectability  thresholds 
based  on  noise  level  and  station  magnitude  corrections  was 
only  partly  successful.  The  results  were  generally  consistent 
with  previous  estimates,  but  with  occasional  deviations  of 
several  tenths  of  an  m^  unit.  The  simple  method  of  estimating 
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station  detecting  capability  by  counting  the  number  of 
reported  events  was  found  to  give  a good  first  order  estimate 
of  relative  50  per  cent  thresholds  between  stations. 

For  each  station,  there  are  significant  regional  variations 
in  detection  thresholds  that  cannot  be  attributed  to  dif- 
ferences in  epicentral  distances.  Of  particular  interest  is 
the  observation  that  the  50  per  cent  threshold  on  the  average 
is  0.4  mfa  units  higher  for  Western  North  America  than  for  the 
Japan-Kuriles-Kamchatka  arc,  although  the  epicentral  distances 
are  similar  for  these  two  regions.  This  points  to  a consider- 
ably stronger  attenuation  of  high  frequencies  in  the  mantle 
for  the  paths  corresponding  to  the  former  region. 

An  indirect  argument  was  used  to  compare  the  performance  of 
tae  Fennoscandian  stations  to  the  average  ISC  station  for 
ore  selected  region  (Japan-Kuriles-Kamchatka) . It  appears 
that  the  Fennoscandian  stations  are  generally  better  than 
average  for  this  region,  with  a mean  improvement  in  50  per 
cent  thresholds  of  about  0.4  mb  units.  A more  complete  evalua- 
tion of  the  relative  standings  of  individual  stations  in  the 
world-wide  network  would  require  a much  more  extensive 
study,  but  could  possibly  be  obtained  using  the  general 
techniques  described  previously  in  this  paper. 

S.  Pirhonen  (Helsinki) 

F.  Ringdal 

K.A.  Berteussen 
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VII. 7 


On  the  M_-m^  Relationship  of  Earthquakes 

The  purpose  of  this  study  is  to  estimate  the  Mg-mb  relationship 
using  data  from  VLPE,  ALP A,  NORSAR  and  PDE,  and  assess  the 
impact  of  network  magnitude  bias  effects  (Ringdal,  1976)  on 
the  results.  We  assume  that  there  exists  a linear  functional 
relationship  between  Mg  and  mfa  of  shallow  earthquakes  of 
the  form 

Ms  = a * mb  + 0 + e (1) 

where  a and  6 are  unknown  coefficients  and  e is  a random 
variable  that  is  normally  distributed  with  zero  mean  and  an 
unknown  variance.  Generally,  the  formula  (1)  can  only  be 
used  as  an  approximation  over  limited  magnitude  ranges;  for 
example,  Gutenberg  and  Richter  (1956)  estimate  a=1.6  using 
only  very  large  events,  while  Evernden  (1975)  finds  a value 
of  ot=1.0  below  magnitude  mb=5.0  down  to  at  least  mb*3.0. 

This  study  is  restricted  to  the  magnitude  ranqe  of  most  inter- 
est for  current  seismic  discrimination  studies,  i.e.,  aporoxi- 
mately  mb~4.0  to  6.0.  The  data  base  consists  of  52  randomly 
selected  shallow  Eurasian  earthquakes  (Turnbull  et  al,  1975) . 
Magnitudes  of  these  events  have  been  available  from  PDE  and 
NORSAR  (mb),  and  from  VLPE,  ALPA  and  NORSAR  (Mg) . In  addition, 
we  used  the  maximum- likelihood  technique  of  Rinqdal  (1976) 
to  modify  the  PDE  mb  estimates,*  we  denote  these  modified 
estimates  by  PDE(m-l)  i»b  values.  The  assumptions  required 
to  obtain  the  PDE(m-l)  estimates  are  described  in  detail  by 
Turnbull  et  al  (1975). 

In  order  to  examine  the  variations  of  the  Mg-n^  slope  as 

a function  of  estimation  techniques  of  M and  m.,  a total 

s D 

of  nine  cases  were  run  based  on  the  given  event  population. 

In  each  case,  the  M values  estimated  by  either  VLPE 

s 

(averaging),  ALPA  or  NORSAR  \re re  combined  with  the  m.  values 

o 
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of  NORSAR,  PDE  and  PDE(m-l).  One  case  is  shown  in  Fiq.  vil  7 l 
while  complete  results  are  summarized  in  Table  VI I. 7.1.  The 
following  points  are  noteworthy: 

1.  Four  of  the  runs  produce  a virtually  identical  slope 
(a^i. 37).  These  are  precisely  thosi*  four  runs  that  combine 
"consistent"  Mg  and  mfc  values,  i.e.,  values  free  of 
network  bias  effects.  (NORSAR  or  ALK'v  Mg  versus  NORSAR 

or  PDE(m-l)  mfa) . 

2.  A consistently  high  value  of  the  slope  (1.66  or  1.64)  is 

found  when  PDE  m.  is  plotted  against  a consistent  M . 

s 

3.  A consistently  low  value  of  a (1.23  or  1.24)  results 
when  VLPE  M is  plotted  against  a consistent  m.  . 

D 

4.  When  PDE  mfa  is  plotted  against  VLPE  Mg,  a is  again  high, 
showing  that  the  network  bias  effects  in  PDE  magnitudes 
dominate  those  of  VLPE. 

Hence,  the  behavior  of  the  computed  slope  aqrees  well  with  what 
could  be  expected  from  network  bias  considerations  (network 
magnitude  are  expected  to  produce  a bias  that  is  largest  for 
small  events) . It  appears  that  the  most  accurate  linear 
functional  relationship  between  Mg  and  mfa  for  the  given 
data  set  (ranging  in  m^  values  from  about  4.0  to  6.0)  has 
a slope  of  approximately  1.4. 

Considering  more  closely  the  four  cases  of  consistent  esti- 
mates, it  is  interesting  to  note  that  the  value  of  the 
orthogonal  standard  deviation  o is  lower  when  using  PDE(n-l) 
mb  versus  either  ALPA  or  NORSAR  Mg  compared  to  when  NORSAR  mfa 
is  used  (ow0.26  vs  o« 0.31).  It  would  be  interestino  to  compare 
PDE(m-l)  mfa  to  the  VLPE-ALPA-NORSAR  combined  network  with 
the  network  bias  reduced  by  maximum  likelihood  processing,* 
however,  we  have  not  been  able  to  do  this,  mostly  because  of 
the  lack  of  reliable  VLPE  data  for  some  stations. 

F.  Ringdal 

L.S.  Turnbull  (Texas  Instruments, 

USA) 
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Table  VII. 7.1 


Estimated  relationship  Ms_m^  for  various  combinations  of  Mg  and  mfa 

measurement  procedures. 


“b 

a 

s 

No. 

of 

Points* 

NP 

Slope 

a 

Intercept 

0 

0 

Orthogonal 

VLPE 

52 

1.54 

-3.58 

■H' 

PDE 

ALPA 

41 

1 . 6l 

-4.33 

NORSAR 

35 

1.64 

-3.95 

0.276 

PDE 

VLPE 

52 

-1.73 

K 

Maximum 

ALPA 

41 

-2.57 

mESm,: 

Likelihood 

NORSAR 

35 

1.37 



-2.25 

0.259 

VLPE 

52 

m 

fl 

WSm 

NORSAR 

ALPA 

41 

NORSAR 

35 

BSP 

B 

The  missing  data  points  (NP<52)  are  due  to  lack  of  available 
data  for  NORSAR  or  ALPA  M for  some  events,  and  not  due  to 
nondetection  at  these  stations. 


An  Improved  Discriminant  using  Spectral  Estimates  for 
Surface  Waves 

In  Bungum  and  Tj0stheim  (1976)  an  improvement  of  the  m^rM 
discriminant  was  obtained  by  combining  the  mb:Mg  data  with 
some  parameters  resulting  from  modelling  the  short  period 
data,  noise  as  well  as  main  signal  and  coda,  as  3rd  order 
autoregressive  time  series.  The  new  two-dimensional  Z1:Z2 
discriminant  was  defined  as 

Z1  = mb  - 0.4  ax (S) 

(1) 

Z2  = Ms  + 0.4  (ax (C)  - (N) ) 

where  a^S)  is  the  first  order  autoregressive  coefficient 
from  the  main  signal  and  C and  N denote  coda  and  preceding 
noise  respectively. 


We  have  now  investigated  the  effect  of  replacing  the  surface 
wave  magnitude  M by  a long  period  power  spectrum  estimate. 

More  precisely , we  adopted  ^20'  the  value  of  the  estimated 
spectral  density  (using  maximum  entropy  spectral  estimation 
over  suitable  group  velocity  windows)  at  a period  of  20  sec 
as  an  estimate  of  the  incoming  energy.  Using  correction  formulae 
for  ordinary  Mg  measurements,  we  constructed  the  distance 
corrected  energy  estimates 

E20  = 7 lo^10l^20/7l2)  + i‘66  1o9i0(A) 

for  distances  A larger  than  25°  and 

E20  = 7 log10(*20/,,2)  + 1o9io(A)  + °*92 

for  distances  A less  than  25°.  Vhese  quantities  were  computed 
for  a Eurasian  data  set  consisting  of  43  earthquakes  and  46 
explosions.  This  data  set  is  identical  to  the  main  data  set 


of  Bungum  and  Tj0stheim  (1976)  with  the  exception  of  two 
earthquakes  for  which  no  digital  long  period  data  could  be 
found.  The  parameters  an^  E2Q  were  computed  not  only 
for  Rayleigh  (RV)  vertical  component  waves  but  also  for 
Rayleigh  (RH)  radial  component  waves  and  Love  (L)  transverse 
component  waves.  Furthermore,  a number  of  different  group 
velocity  data  windows  were  used.  The  information  of  interest 
for  discrimination  purposes  is  summarized  in  Table  VII. 8.1, 
where  mean  values  have  been  computed  for  three  of  the  windows 
studied.  On  the  basis  of  the  information  in  this  table  we 
constructed  the  quantities 

Y2  = (Love)  + (E^o*  (Love>  - E 2 q ^ ( Love ) ) (1) 

where  the  superscripts  refer  to  the  windows  used.  Replacing 
Mg  with  Y2  we  arrive  at  the  X1:X2  discriminant  given  as 


XI  = mfa  - 0.4  a, 

X2  = Y2  + 0.4(a, (C)  - ax (N) ) 


(2) 


Our  results  indicate  that  the  discriminant  (2)  works  sig- 
nificantly better  than  the  discriminant  in  Bungum  and  Tj0stheim 
(1976).  In  Figs.  VII. 8.1  and  VII. 8. 2 the  discriminant  (2)  is 
compared  with  the  mfa:Mg  discriminant.  The  mb:Mg  diagram  has 
two  explosions  within  the  earthquake  population  and  the 
separation  between  earthquakes  and  explosions  is  not  particu- 
larly good.  In  the  XliX2  diagram  these  two  events  are  now 
moved  well  into  the  explosion  population  and  the  separation 
between  the  two  populations  is  substantially  improved. 


Dag  Tj0stheim 
Eystein  Husebye 
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Table  VII. fl. 1 

The  mean  values  of  the  energy  estimates  E«  and  of 
surface  wave  magnitude  Ms  as  obtained  by  averaging  over  the 
46  explosions  and  43  earthquakes  respectively. 


1 

2 

3 

Group  Velocity  Windows  (km/s) 

8. 0-2. 2 

3. 3-2. 7 

E2Q  Rayleigh  vertical,  explosions 

3.88 

3.89 

— 1 

3.61 

^20  Raylei9h  vertical,  earthquakes 

4.40 

4.56 

4.05 

E2Q  Rayleigh  horizontal, 

explosions 

3.85 

3.81 

3.54 

E2q  Rayleigh  horizontal, 

earthquakes 

4.40 

4.47 

3.99 

Group  velocity  windows  (km/s) 

8. 0-2. 2 

3.7-3. 1 

4 .6-3.8 

E20  Love,  explosions 

3.77 

3.76 

3.51 

E2o  Love,  earthquakes 

4.55 

4.69 

4.00 

M , explosions 
s 

3 72 

Mg,  earthquakes 

4.39 

5UR‘cHLt  UP/E  MAGNITUDE  MS 


2 3 H 5 B 

BODY  NAVE  MAGNITUDE  MB 


Fig.  VII. 8.1 


mb:Ms  diagrajn  for  the  46  explosions  and  43  earthquakes. 

PDE  nl  and  NOR ShR  M values  (rounded  to  the  first  decimal) 
have  Been  used.  The  regression  lines  have  been  fitted  using 
the  same  maximum  likelihood  procedure  as  in  Bungum  and 
Tj0stheim  (1976).  y 


3GEIMIMRNT  PARAMETER  X2 
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DISCRIMINANT  PARAMETER  XI 


Fig.  VII. 8. 2 ^a9rara  ^or  t^ie  46  explosions  and  43  earthquakes. 


VII. 9 The  Seismicity  of  Fennoscandia 

A critical  parameter  in  seismicity  studies  is  the  space-time 
relationship  between  the  largest  occurring  earthquakes 
within  the  region  under  investigation.  An  interesting  feature 
here  is  that  most  of  the  largest  shocks  (M>6)  are  restricted 
to  the  period  1863-1913,  while  3 other  strong  earthquakes 
took  place  between  1819  *?nd  1836.  There  is  no  indication  of 
spatial  migration  in  the  distribution  of  the  largest  earth- 
quakes and  neither  do  we  find  any  clear  regularity  in  occur- 
rence of  large  events  in  specific  areas.  In  other  words,  the 
very  largest  earthquakes  observed  seem  to  be  isolated  phenomena 
This  points  towards  exceptionally  long  recurrence  intervals, 
i.e.,  a large  time  lag  between  two  lazge  events  occurring 
within  a tectonically  uniform  area.  As  it  is  considered  un- 
likely that  any  major  earthquake  occurring  after  1600  are  left 
unreported,  the  recurrence  interval  is  probably  larger  than 
350  years  for  earthquakes  with  magnitude  M larger  than  6. 


The  seismicity  pattern  of  Fennoscandia  is  somewhat  diffuse, 
which  is  normal  for  intraplate  earthquake  occurrence  (Sbar 
and  Sykes,  1973) . We  found  it  convenient  to  subdivide  the 
Fennoscandian  earthquakes  in  3 primary  zones:  namely,  the 
western  Norv’ay  zone,  the  Telemark-Vanern  zone  and  the  Bethnian 
zone,  plus  a weaker  one  called  the  Lappland  zone  (Fig.  VII. 8.1) 
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The  western  Norway  seismicity  belt  is  in  northern  Norway 
confined  to  the  coastal  area  while  the  epicenters  become 
more  dispersed  south  of  64°N.  The  west  coast  area  of  Norway 
is  relatively  prominent  seismically  and  also  geophysically 
(Husebye  et  al,  1975).  A characteristic  feature  here  is  that 
this  zone  is  within  the  Caledonides  and  its  strike  direction 
is  parallel  to  the  folding  axis)  this  can  be  taken  as  evidence 
for  the  importance  of  remnant  or  loc  ed-in  stresses  from  the 
above  mountain  folding  period.  Another  feature  is  that  the 
epicenters  are  mainly  confined  to  the  coastal  areas  where 
the  relief  is  very  pronounced,  which  in  turn  may  indicate 
additional  loading  stresses  and/or  zones  of  weakness. 

The  Telemark— Vanern  seismicity  zone  which  represents  a 
geographical  envelope  of  earthquake  epicenters  in  this  area 
is  typified  by  graben  structures  and  the  most  prominent  one 
is  the  Oslo  Graben.  A slight  clustering  cf  epicenters  can  be 
found  around  the  mentioned  graben  structures,  pointing  towards 
a causal  connection  between  these  structures  and  earthquake 
occurrence,  possibly  through  release  of  locked-in  stresses. 

An  alternative  explanation  is,  however,  that  the  orientation 
of  the  Telemark-Vanern  belt  is  in  fact  tectonically  significant. 

The  Bothnian  seismicity  belt,  parallel  to  the  Caledonian  folding 
axis,  goes  from  Lake  Vattern  to  the  northern  end  of  the  Gulf 
of  Bothnia.  The  latter  area  has  the  most  pronounced  earthquake 
activity,  and  there  may  be  a correlation  between  the  relatively 
strong  glacial  uplift  and  the  earthquake  occurrence  in  this 
particular  area  (BSth,  1953).  ALso,  this  zone  is  characterized 
by  faulting. 


The  Lappland  seismicity  zone  is  only  weakly  defined  by  the  data 
available  to  us,  but  a series  of  earthquakes  in  1973/74  was 

entirely  within  this  zone.  As  this  area  always  has  been 
thinly  settled,  the  macroseismic  information  available  would 
necessarily  be  scarce. 
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The  above  seismic  zones  account  for  most  of  the  seismic  activity 
in  Fennoscandia  during  the  last  five  hundred  years.  The  seismic 
activity  is  typical  for  intraplate  earthquake  occurrence,  by 
neither  exhibiting  a too  clear  spatial  zoning  nor  an  obvious 
correlation  with  geological  and  geophysical  information  per- 
tinent to  the  area.  The  likely  reason  for  this  is  that  Fenno- 
scandia  has  been  through  several  tectonic  cycles,  which  in  turn 
is  reflected  in  the  present  complex  stress  distribution.  Con- 
sequently, an  improved  understanding  of  the  on-going  seismic 
activity  here  and  at  the  same  time  a better  assessment  of 
dominant  stress  sources  requires  more  in- situ  stress  measure- 
ments. Equally  important  would  be  focal  mechanism  solutions 
for  earthquakes  occurring  within  Fennoscandia. 

E.S.  Husebye 

H.  Bungum 
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Fennoscandian  earthquakes  for  the  time  period  1497-1073  and 
with  a magnitude  M greater  than  4.5.  Ar  outline  is  also 
given  of  the  4 seismicity  zones  defined  in  this  paper.  Note 
that  the  definition  of  the  zones  is  based  on  much  more 
seismicity  data  than  is  actually  displayed  in  this  figure. 
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VII. 10  Precise  Monitoring  of  Seismic  Velocities 

The  physical  setup  for  this  experiment  was  explained  in  the 
previous  report  in  this  series  (NORSAR  Scientific  Report  No. 
-’”74/75).  ■‘■t  involves  a power  station  as  a source  vibrating 
at  frequencies  around  2.778  Hz,  5 seismometers  (Ch  1-4 
and  6 at  n ’array  14C) , at  distances  from  4.7  to  13.7  km 
from  the  source,  and  a 6th  data  channel  (Ch  5)  where  the 
network  power  of  50  Hz  is  input  after  being  divided  18  times 
down  to  2.70  Hz  and  properly  reduced  in  voltage  (to  +4.6  V). 

An  extersive  analysis  has  now  been  completed  involving  about 
120  hours  (5  days)  of  data.  We  have  Fourier-analyzed  altogether 
425  blocks  each  of  10  000  samples  of  10  Hz  data  (16  2/3 
minutes  per  block),  giving  a frequency  resolution  of  0.001 
Hz.  For  each  block,  the  reference  channel  (Ch.  5)  was  first 
analyzed  so  as  to  find  out  what  frequencies  to  analyze, 
which  is  important  because  the  network  power  frequency 
actually  fluctuates  as  much  as  +0.3%  around  the  nominal 
frequency,  with  a correlation  time  of  sometimes  only  a few 
seconds.  This  corresponds  to  frequencies  in  the  ranqe  2.771 
to  2.785  Hz,  where  our  resolution  gave  15  independent  esti- 
mates. Tor  each  block,  we  then  analyzed  any  frequency  in 
that  range  at  which  the  reference  channel  had  stayed  for  at 
least  7%  of  the  time,  usually  givinq  between  3 and  6 fre- 
quencies. Having  established  the  frequencies  to  analyze  in 
this  way,  the  initial  analysis  involved  estimating  the 
Fouvier  transform  for  each  of  the  5 seismic  data  channels 
for  each  of  the  selected  frequencies.  These  basic  Fourier 
Transform  results  were  then  stacked  on  a digital  tape,  ready 
for  subsequent  analysis. 

The  way  to  monitor  seismic  velocities  from  these  data  is 
through  the  analysis  of  spectral  phase  differences.  We  will 
here  only  discuss  3 of  the  many  possible  combinations,  namely, 
between  the  source  (Ch  5)  and  Ch  4 (distance  13.7  km),  between 
the  source  and  Ch  1 (distance  4.7  km),  and  between  Ch  1 and 
Ch  2.  The  latter  two  channels  are  the  ones  closest  to  the 
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source  and  therefore  with  the  best  signal-to-noise  ratio, 
they  have  also  practically  the  same  distance  to  the  source 
and  are  consequently  essential  as  a control  combination. 

The  phase  estimates  in  Fig.  VII. 10.1  are  obtained  by  summing 
the  power  spectral  estimates  over  all  the  analyzed  blocks  of 
data,  after  some  reduction.  For  each  of  the  frequencies  in 
each  of  the  blocks  ve  analyzed  the  power  in  the  reference 
channel.  It  has  been  found  that  the  quality  of  the  phase 
estimates  is  proportional  to  this  power  level,  and  conse- 
quently a threshold  was  determined  below  which  the  data  was 
rejected.  In  this  way  we  improved  the  quality  of  the  phase 
estimates  by  rejecting  about  40%  of  the  original  data,*  in 
Fig.  VII. 10.1  it  is  shown  that  there  remains  between  1 and 
146  blocks  for  each  of  the  frequencies. 


An  important  by-product  of  the  phase  estimates  in  Fig.  VII. 10.1 
is  that  they  actually  give,  albeit  with  considerable  uncer- 
tainty, an  estimate  of  the  group  velocity,  either  between 
two  seismometers  or  between  any  seismometer  and  the  source. 

The  slope  d<J>/du)  has  been  estimated  using  a weighted  least 
squares  estimation  procedure,  using  as  weights  the  number 
of  blocks  for  each  frequency  (since  that  number  is  directly 
proportional  to  the  stability).  Once  the  slope  is  estimated, 
the  group  velocity  is  simply  the  distance  divided  by  the  slope. 
The  group  velocity  results  are  not  really  as  consistent  as 
could  be  the  impression  from  Fig.  VII. 10.1  (v  = 3.9+0.b  and 
3. 8+0. 4 km/s),  for  the  combinations  of  channels  with  accept- 
able signal-to-noise  ratios  the  estimates  actually  vary 
between  3.0  and  5.0  km/s.  However,  using  the  group  velocities 
arrived  at  in  this  way,  and  assuming  the  same  values  for  the 
phase  velocities,  it  is  possible  to  estimate  the  total 
f ravel  time  and  thereby  the  relative  size  of  observed  phase 
differences.  It  is  fourd  in  this  way  that  a precision  of  +10-^ 
for  combination  1-5  corresponds  to  +1.2°,  and  for  4-5  the  value 
is  +3.6°,  corresponding  to  the  error  bars  in  Fig.  VII. 10.1. 
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An  obvious  way  to  further  improve  the  stability  of  the  phase 
estimates  is  to  use,  instead  of  the  values  for  individual 
frequencies,  the  intercept  of  the  regression  line  with  the 
center  frequency.  (The  error  bar  in  Fig.  VII. 10.1  has  been 
placed  around  that  value.)  That  is  the  technique  used  in 
Fig  VII. 10.2,  where  each  phase  estimate  is  the  average  over 
200  minutes  of  data.  On  the  average,  there  are  around  30 
individual  spectral  estimates  (12  blocks  and  2.5  frequencies 
per  block)  behind  each  data  point  in  Fig.  VII. 10. 2.  The 
figure  shows  that  there  is  no  significant  time  variation 
in  the  data,  and  the  standard  deviation  for  combinations 
1-5  and  4-5  are  1.5°  and  4.0°,  respectively,  which  is  quite 
cloie  to  the  values  given  above  for  a precision  of  10~\ 

The  fact  that  we  get  the  same  precision  for  the  two  combina- 
tions means  that  the  loss  in  signal-to-noise  ratio  for  the 
greatest  distance  is  compensated  Ly  the  increased  distance 
(ay  a factor  of  3).  Assuming  now  that  a precision  of  10*^ 
is  obtained  for  each  of  the  points  in  Fig.  VI 1. 10. 2,  we  should 
by  summing  all  the  data,  obtain  a precision  of  about 

10  //I&*2«10  , which  refers  to  the  intercept  points  In 

Fig.  VII. 10. I. 

This  report  is  only  preliminary,  and  work  continues  along 
the  following  lines:  1)  Predict  the  possible  influence  from 
tidal  stress  variations  (the  5 days  analyzed  here  are  selected 
so  as  to  cover  a time  period  when  those  variations  are  sup- 
posed to  be  at  a maximum),  2)  Further  improve  (by  proper 
data  reduction)  and  more  accurately  assess  (from  spectral 
variances)  the  precision  of  the  already  acquired  data  base, 

3)  As  a continuation  of  the  project,  consider  using  higher 
harmonics  of  the  2.778  Hz  data,  which  have  quite  good  signal- 
to-noise  ratios,  and  where  tiv«  phase  differences  due  to  tidal 
effects  should  be  larger. 


H . Bungum 

T.  Risbo  (Copenhagen) 

E.  Hjortenberg  (Copenhagen) 
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